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We show that the physical system consisting of trapped ions interacting with lasers may undergo 
a rich variety of quantum phase transitions. By changing the laser intensities and polarizations 
the dynamics of the internal states of the ions can be controlled, in such a way that an Ising or 
Heisenberg-like interaction is induced between effective spins. Our scheme allows us to build an 
analogue quantum simulator of spin systems with trapped ions, and observe and analyze quantum 
phase transitions with unprecedent opportunities for the measurement and manipulation of spins. 



o 
o 

OS 



> 
O 



o - 
O ' 

^ : 



Introduction. Quantum spin models are a paradigm 
for the study of many-body effects, and show intriguing 
phenomena like the existence of quantum phase transi- 
tions at certain values of the parameters that govern the 
spin Hamiltonian 1]. In this work we show that quantum 
phase transitions can be induced by lasers in ion traps, 
where the internal states of the ions play the role of ef- 
fective spins. In our scheme an Ising, XY, or XYZ spin- 
spin interaction is transmitted by collective vibrational 
modes and can be switched and tuned by the lasers and 
by the choice of trapping conditions. The understand- 
ing of quantum spin models has been hindered by the 
fact that they are very hard to be numerically simulated. 
This work presents a way to circumvent this difficulty 
by using a system of trapped ions as an analog quantum 
simulator for magnetic systems 0, l^- Other schemes 
for the simulation of quantum problems [4|, IM lal rely on 
the use of a quantum computer 0, 0, 1^ lift fill, or on 
the stroboscopic change of quantum parameters |l2l Il3j | . 
Our proposal is not based on quantum gates, so that the 
requirements for its implementation are much less strin- 
gent, and it is not exposed to the errors accumulated in 
stroboscopic methods. 

The main advantage of ions is that they can be trapped 
and cooled very efficiently, and therefore they can be 
stored at fixed positions in space |llj|. Furthermore, their 
internal states can be precisely manipulated using lasers, 
and measured with basically 100% efficiency [l5|- Our 
proposal, thus, would make feasible to study quantum 
phase transitions in this particular experimental system, 
where one is able to tune the range and strength of the 
interaction, and manipulate and measure single effective 
spins, something that is not possible in solid-state mag- 
netic systems. 

Spin-spin interaction. Our scheme works with a set 
of N ions trapped by electric and/or magnetic forces. 
For example, they could be ions in a linear Paul trap, in 
microtraps in 1 or 2 dimensions, or in a Penning trap. 
In order to simplify the discussion that follows, we fo- 
cus here on the one dimensional case: a chain of trapped 
ions. Let us assign the z (a;, y) vectors to the axial (ra- 
dial) directions relative to the chain. We will assume that 
the ions have two internal ground hyperfine levels, which 
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FIG. 1: (color online) Possible laser configuration such that 
the ions experience a force along the a = x, y, and z di- 
rections, if they are in the state | |)a (Fig (a)). Each pair of 
copropagating lasers independently excite a Raman transition 
IT) ~^ I i) with different detunings from the virtual excited 
state |e) (such that they do not interfere). The relative phase 
of each pair of copropagating lasers along the direction a is 
chosen such that the state | J,)a is dark. Thus, they effectively 
interact only with the level \1)a [see Fig (b)]. Counterpropa- 
gating lasers form a standing wave pattern which pushes the 
ions along the direction a only if they are in state | |)c,. Note 
that along the z direction only one laser is needed. Of course, 
this direction does not need to coincide with the trap axis for 
the case of a linear ion trap. 



play the role of components of an effective spin 1/2 (our 
scheme can be easily extended to simulate larger spins). 
We denote by ct", a = a;, y, z, the Pauli operators, and by 
I J,)q and 1 1)a the corresponding eigenstates. All the ions 
are driven by the same off resonant laser beams prop- 
agating along the three spatial directions (Fig. QJ. In 
the most general configuration, we will assume that the 
lasers propagating along a push the ions in that direction 
provided they are in states | t)a- This can be achieved 
by choosing the relative phases of the lasers so that the 
state \l)a is dark with respect to the lasers propagating 
along the direction a (see, e.g., Fig.^. 

As a result of Coulomb repulsion, vibrational de- 
grees of freedom are collective modes (phonons): Hy 



I, where a]j „ {liu)a,n): are the opera- 
tor (energy) of the n mode in the a direction. Our Hamil- 
tonian includes also the force produced by the lasers 
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(Hf), and the action of effective magnetic fields {Hm)'- 
H = Hy + Hf + Hm ■ 



the equihbrium positions along the chain, in such a 
way that the Coulomb repulsion is balanced with the 
trapping forces. The second derivatives of the Coulomb 



ff^ — _2 "y^Faqa i| t) (T \a ij H „i ^ B"" a"" . (1) energy with respect to these displacements determine the 



elastic constants of the chain 



The fields B" can be simulated by lasers acting on the 
internal transition of the ions only, qa^i are the three 
spacial coordinates of each ion, and i refers to the site 
in the chain. The local coordinates can be expressed in 
terms of the collective modes: 



Hf 



(2) 

A suitable method to get the promised spin Hamilto- 
nian is to apply the following canonical transformation: 



(3) 



In the new basis the Hamiltonian includes an effective 
spin-spin interaction: 

e-^i?e^ = iJ, + i ^ J^af + ^ B'Vf + He, (4) 



where 



The effective magnetic fields receive a contribution from 
the pushing forces -B'" — B°' + F^/(mtj^), where Wq are 
the trapping frequencies in each direction. Note that the 
extra term in B'" does not depend on the site of the ion. 
He describes the residual coupling between the effective 
spins and the collective vibrational modes. To lowest 
order in 77"^ , it is given by the following expression: 



H. 



(6) 



In the transformed basis, Hf disappears, and we get an 
anisotropic Heisenberg (XYZ) interaction between effec- 
tive spins. He is a perturbation that limits the accuracy 
of the simulation, and can be neglected if the adimen- 
sional coefficients 77f„ are small enough. We will discuss 
later this approximation and consider now the properties 
of the effective spin-spin interaction. 

J^j depends on the characteristics of the vibrational 
modes in the corresponding a direction. The ions occupy 
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where Cx^y = Cz — —2. The unitary matri- 
ces in (0) diagonalize the vibrational hamiltonian: 
Mfj^K.f.jMj.^ = ujl „^Sn,m, and Eq. (01 implies that we 
can express the effective interaction in the following way: 
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If we assume a constant distance between ions, do, then 
we can summarize the properties of the vibrational modes 



by means of the parameters Pq 



Ca e- 



'■/muj^dl. We 



distinguish two cases: (i) Stiff modes: 1 im- 

plies that the Coulomb interaction can be considered as 
a perturbation to the trapping potential, so that the 
inverse of the elasticity matrix /C" can be calculated 
to first order in the non-diagonal terms in Eq. JT)): 
J"_^j — CaF'^e^ / {muj1^Y\z^ ~ -^il^- ^^^^ '^^^^ S*^* 
a spin-spin interaction with a dipolar decay law, that 
is antiferromagnetic (ferromagnetic), if it is transmitted 
by the radial (axial) vibrational modes, (ii) Soft modes: 
/3q ^ 1 implies that the Coulomb interaction is impor- 
tant against the trapping of the ions in the a direction, 
a situation that is discussed below in connection with 
concrete trapping conditions. 

Experimental realizations. We discuss now the experi- 
mental setups that can be used to implement the simu- 
lation of the quantum spin hamiltonians exposed above. 

(i) Arrays of ion microtraps. In this case ions are con- 
fined individually by microtraps, in such a way that they 
form a regular array in 1 or 2 dimensions 'iS'. A ID ar- 
ray of microtraps can be realized, for example, by using a 
linear Paul trap, and an additional off-resonant standing- 
wave that creates a periodic confining potential in the ax- 
ial direction. In the ID case, forces in the three spacial 
directions would simulate an effective XYZ interaction, 
antiferromagnetic, or ferromagnetic, for the (cc, y), or 
z spin components, respectively. If the radial trapping 
frequencies large enough, the stiff limit could be 

reached, and the spin-spin interaction would have a dipo- 
lar decay. By relaxation of the radial trapping frequencies 
the spin-spin interaction could be controlled and, in par- 
ticular, the interaction between second neighbors could 
be suppressed, similarly to the case of linear Paul traps 
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FIG. 2: (color online) Effective spin-spin interactions in a 
Coulomb chain. Left: long-range ferromagnetic interaction 
transmitted by the longitudinal modes (A'^ = 50 ions). Right: 
short-range antiferromagnetic interaction transmitted by the 
radial modes (N = 100). For /3x,y = 0.01 (lower curves), 
the effective interaction shows a dipolar decay: Ji.j oc |i — 
, while for Px,y ~ 0.1 (upper curves), the second-nearest 
neighbor interaction is suppressed. 

(see below). The effective XYZ interaction, with the ad- 
dition of effective magnetic fields would allow to study a 
rich variety of quantum phase transitions. 

(ii) Linear Paul traps. Paul traps are currently used to 
cool and arrange ions in ID Coulomb crystals, in which 
the distance between ions is not constant. We can, how- 
ever, define an averaged lattice constant do, in order to 
understand the qualitative properties of the vibrational 
modes. We distinguish two cases, depending on the ori- 
entation of the pushing forces relative to the chain: (a) 
Axial force. The equilibrium position of the ions are such 
that (3z ^ 1 01 . The center of mass mode has the main 
contribution in Eq. ©, so that we get a long-range fer- 
romagnetic interaction. Our numerical calculation for a 
chain of = 50 ions (Fig. [ll(left)), shows that a force in 
the z direction results in an interaction, Jij, with a range 
comparable to the size of the chain, (b ) Radial force. The 
characteristics of the effective spin-spin interaction dras- 
tically change when lasers push the ions in the radial (x, 
y) directions. The stability of the chain against zig-zag 
deformation implies that f3x,y < 1 |l7l |. If the trapping 
frequencies lOx, Wy are large enough, the limit in which 
the spin-spin interaction shows a dipolar decay can be 
easily reached. For P^.y ^ 0.05, the spin-spin interaction 
departs from the dipolar case, and the main contribution 
comes now from the collective stretched mode. With 
the proper choice of trapping frequencies, the second- 
neighbor term in the dipole interaction can be suppressed 
(see Fig. |21 (right)). The possible experimental configu- 
rations in a linear Paul trap are the same as those for 
a linear array of ion microtraps, with the difference that 
the modes in the axial direction transmit a long-range 
interaction. 

(Hi) Penning traps. Finally, 2D Coulomb crystals 



could be used for the simulation of 2D quantum spin 
systems. Cooling of ions in a 2D hexagonal lattice has 
been demonstrated in Penning traps 0, . This ex- 
perimental setup has the additional advantage that the 
hexagonal lattice could allow to simulate the effect of 
magnetic frustration. For example, a pushing laser in 
the direction perpendicular to the 2D crystal would in- 
duce an antiferromagnetic Ising interaction. The decay 
of this interaction would also be dipolar, if the 2D vibra- 
tional modes transverse to the crystal plane are in the 
stiff limit. 

In all the three experimental systems discussed above 
the simplest laser configuration would consist of only 
one force in a given direction, so that an Ising interac- 
tion cx Ja^a^ would be simulated. With the addition 
of an effective transverse field oc Ba^, the spin system 
would present a quantum phase transition between a fer- 
romagnetic (or antiferromagnetic) ordered ground state 
(J B^), and a disordered ground state {B'^ :S> J), 
at a value J « B^ Q. For example, one could proceed 
as follows: (1) Prepare the state | i ... i). (2) Switch 
on adiabatically the magnetic field B^. (3) Switch on 
adiabatically the effective interaction up to a given value 
J. (4) Measure the global fluorescence. This sequence 
could be repeated with several values of J, starting with 
J <^ B^ , and increasing the value of the spin-spin inter- 
action in each experiment. The phase transition would be 
evidenced in the global fluorescence, that would show the 
population of the internal states of the ions, and, thus, 
the emergence of ferromagnetic (or antiferromagnetic) or- 
der in the effective spins. Furthermore, the possibility of 
addressing single ions could allow us to manipulate the 
many-body ground state, create excitations (spin- waves) , 
or study the entanglement between effective spins. The 
non-equilibrium quantum dynamics of the Ising model 
could also be studied, for example, by switching abruptly 
the interaction on. 

Validity of the effective spin Hamiltonian. Assume 
that the ions are initially in the internal state {ipi), that 
would evolve to \^f) = exp{~iHst/h)\'iJji) , under the ac- 
tion of the simulated spin Hamiltonian, Hs. The fidelity 
of the simulation is given by the probability of finding 
\ipf) after the real evolution: 

Ti^:,) = {e-'"'/^ ® pph e'"'/^)i,f = (9) 

where pi = |i/;i)('0i|,and the vibrational modes are ini- 
tially in a mixed thermal state (pph). Hq = Hy + Hs, 
is the Hamiltonian without residual spin-phonon cou- 
pling. In order to simplify the discussion, we define 
rj ~ F y'h/Qmuj / huj {a is omitted), as the parameter that 
characterizes the displacement of the vibrational modes, 
and estimate the error, E = 1 — J-, in a, series in rj. If we 
use only one force in a given direction, then He = 0, and 
the error has its origin in the canonical transformation. 
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FIG. 3; (color online) Simulation of the XY model with two 
ions. We have considered the case ujx = IAlu^, and several 
values for ujy. Left: Error of the simulation as a function of the 
mean phonon number n, in a symmetric linear ion trap (uox = 
LOy) (thin line), and in an asymmetric trap, uJy — O.Ybu;^ (thick 
line), with displacements ri = 0.05 (continuous line), and ry = 
0.02 (dashed line). Right: Error of the quantum simulation 
as a function of the displacement of the phonons (77) for the 
asymmetric trap, and n = 0.25. The oscillations are due to 
the non-resonant contributions from each vibrational mode, 
and the averaged error satisfies the relation E (x rj^ (straight 
line). 



A straightforward expansion of in S, allows to esti- 
mate a correction E cx r/^(l -I- 2n), where h is the mean 
phonon number. We estimate now the contribution from 
He- If two or three forces act on the ions, then He gives 
an additional contribution that vanishes at zero temper- 
ature, E (X n. We can avoid the need to cool the ions by 
imposing different trapping frequencies in each direction, 
so that only non-resonant terms appear in Eq. In 
this case, E (x J"^ /{lOx — uiy)'^ — 0(77^), and the main 
contribution to the error is again the one given by the 
canonical transformation. 

We illustrate how the effective spin-spin interaction 
works by presenting a numerical simulation for the case 
of two trapped ions. We consider two lasers that push 
the ions with the same intensity in the radial {x, y) di- 
rections, so that an effective isotropic XY interaction 
(J^ = jy = J) is simulated during a time T — 7r/(8J). 
The evolution of the effective spins under such condi- 
tions corresponds to the application to the initial state 
of a vSWAP gate 5] . The Hamiltonian includes the four 
vibrational modes in the radial directions, and the two 
effective spins that correspond to the internal states of 
the two ions. For the numerical solution, in the range 
or parameters considered here, it is enough to truncate 
the phonon Hilbert space by taking only the three lowest 
states in each mode. As a figure of merit, we have cho- 
sen the averaged fidelity, ^ = f d'il>iT{^i), that describes 
the averaged accuracy of the simulation upon integration 
over all the possible initial states. 

The main results of our numerical calculation are il- 
lustrated in Fig. 13 By introducing a small anisotropy 
{ujx ^ UJy) in the radial trapping frequencies, we get rid 



of the dependence of the error on the mean number of 
phonons, and the requirements for the cooling of the 
trapped ions are not so stringent (Fig. |21(left)). In Fig. 
13 (right) we show the dependence of the error on the dis- 
placement 77, for the case of the asymmetric trap. With 
occupation numbers in the radial modes of fi « 0.25, 
and ujz — 10 MHz, our Eq. |S1 allows us to calculate an 
interaction strength J jYiK, 10 kHz, with error E = 10~^. 

Conclusions and Outlook. We have presented a method 
to simulate Ising, XY, and XYZ interactions between ef- 
fective spins with trapped ions. Our scheme relies on 
the coupling of the internal states with the collective vi- 
brational modes, and would lead to the observation of 
quantum phase transitions in ion traps. With some mod- 
ifications, our proposal could also be used to simulate 
many other systems of great interest in Condensed Mat- 
ter Physics. For example, the action of pushing forces on 
a single ion, could allow to realize a spin-boson model, 
with the bosonic bath being given by the vibrations of 
the ions. Different laser configurations than the one pre- 
sented here, could also allow to simulate other particular 
spin hamiltonians. 
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